Origin of hyperdiffusion in generalized Brownian motion 
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We study a minimal non-Markovian model of superdiffusion which originates from long-range 
velocity correlations within the generalized Langevin equation (GLE) approach. The model al- 
lows for a three-dimensional Markovian embedding. The emergence of a transient hyperdiffusion, 
(Ax 2 (t)) oc t 2+x , with A ~ 1 — 3 is detected in tilted washboard potentials before it ends up 
in a ballistic asymptotic regime. We relate this phenomenon to a transient heating of particles 
Tkin(t) oc i A from the thermal bath temperature T to some maximal kinetic temperature T ma x. This 
hyperdiffusive transient regime ceases when the particles arrive at the maximal kinetic temperature. 

PACS numbers: 05.40.-a, 82.20.Uv, 87.16.Uv 



Anomalous diffusion is met across many different 
branches of physics, from charge transport processes in 
amorphous materials to plasma physics and biophysics, 
and there are many different theories reflecting a variety 
of underlying physical mechanisms [H-0|- One of funda- 
mental approaches is based on the generalized Langevin 
equation (GLE) [8r]23j|. see Eq. (JTJ> below. Remarkably, it 
can be derived from a Hamiltonian dynamics of a particle 
that bilinearly couples to a thermal bath of harmonic os- 
cillators characterized by the bath spectral density J(ui) 
@, IHL El- This model is capable to describe all 
types of anomalous diffusion (Ax 2 (t)) ~ 2D a t a /T(l + a) 
within a unified framework and the index a reflects just 
the low-frequency behavior of th e sp ectral bath density 
J{lo) oc uj a for < a < 2 [H, HJ. Namely, the case 
< a < 1 corresponds to anomalously slow diffusion, or 
subdiffusion; the case a = 1 to normal diffusion, and the 
case 1 < a < 2 to superdiffusion. The ballistic diffusion, 
(Ax 2 (t)) ~ D^t 2 , is attained for all spectral densities 
with a > 2 at low frequencies and the case of J(oj) oc lu 2 
is marginally ballistic (a special case). This finding is 
without a potential, or under a constant force. Then the 
occurrence of hyperdiffusion, (Ax 2 (t)) oc t a with a > 2, 
is not possible if Brownian particles are initially ther- 
malized [25| . Several circumstances are of special inter- 
est. First, GLE diffusion is almost always ergodic, except 
for the ballistic case [13, Qjl US HI] studied also below. 
Similarly, anomalous diffusion based on continuous time 
random walks is weakly non-ergodic [26|. Second, for a 
constant force F, the diffusion coefficient is proportional 
to the bath temperature, i.e. Da oc T and a generalized 
Einstein-Stokes relation holds [12[ • For nonlinear forcing, 
e.g. in tilted washboard potentials, this kind of anoma- 
lous diffusion is not sufficiently investigated and offers 
surprises. In particular, a hyperdiffusive regime occurs 
with a greatly enhanced to a e s ~ 3 — 5 (22|, H3] ■ This 
puzzling nonlinear and nonequilibrium effect is the focus 
of this study. 

Such anomalous diffusion allows for Markovian embed- 
dings of surprisingly small dimensions which suffice nor- 
mally in practice [U [22j. A Markovian embedding is 
natural given that the underlying Hamiltonian dynam- 
ics is Markovian. However, it formally has an infinite 



dimension for a thermal bath considered in the thermo- 
dynamic limit. Surprisingly, the practical embedding di- 
mension using some auxiliary stochastic variables can be 
quite small. 

The purpose of this Letter is to give a physical expla- 
nation of the observed hyperdiffusive anomaly as a tran- 
sient heating of particles with their kinetic temperature 
defined via the velocity variance (Av 2 (t)) oc T^ n rising in 
accordance to a transient power law, T^ n (t) oc f , from 
the bath temperature T to a maximal kinetic temper- 
ature r max which depends on the duration of transient 
period through F, T, and the amplitude of periodic po- 
tential V (T kin = T, when V = 0, or F = 0). It can 
be very large (thousands of T). Such a nonlinear heat- 
ing mechanism in fixed not alternating in time applied 
fields is quite unusual, and, paradoxically, smaller bias 
strengths F yield higher T max values. 

We start from the traditional GLE model in one se- 
lected direction for a particle of mass m in the potential 
V(x) subjected to a linear friction with memory kernel 
?/(i) = (2/n) / °° dujJ(uj) cos{ujt)/uj and random force 
of zero-mean: 

mx+ f n(t - t')x{t')dt' + ^-V{x) = £(t) . (1) 
Jo ox 

The random force is Gaussian and fully characterized 
by its autocorrelation function satisfying the fluctuation- 
dissipation relation 



mat')) = k B Tr 1 {\t-t , \ 



(2) 



which in turn is a consequence of the fluctuation- 
dissipation theorem. 

A necessary condition for the emergence of superdif- 
fusion asymptotically (t — > oo) within the considered 
class of models is zero integral friction |22j, i.e. 
limt-^oo J Q r)(t')dt' = 0. The memory kernel thus must 
be positive at times t' — t, yielding 77(0) > 0, cf. Eq. ((2J, 
and possess a negative part. The simplest model which 
satisfies these two conditions is 1171: 



17(f) = v [25(t) 



(3) 



It will be considered in the following and corresponds to 
a spectral bath density J(uj) which is cubic for oj <C v 
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(typifying e.g. acoustic bulk phonons in solids [l^]) and 
linear for u> 3> v. The corresponding spectral power of 
the noise S(ui) corresponds to the white noise (for lu 3> 
u), i.e. S(lu) = const, with the small-frequency part of 
the spectrum smoothly cut, so that S(oj) oc lu 2 for u> <C v. 

Furthermore, the autocorrelation function (ACF) of 
velocity fluctuations, Av(t) = v(t) — (v(t)), in the ab- 
sence of deterministic force, or under a constant forcing 
obeys for this minimal model 

(Av(t)Av(t')) = v\{^— 
I v + 7 

+_j-exp[-(*, + 7)|*-t'|]} (4) 
v + 7 i 

with 7 = r]/m, provided that the velocities are initially 
thermally distributed with yj (Aw 2 ) — vt — \/k B T ' /m. 
This follows from the Laplace-transformed result for this 
quantity for arbitrary kernels ||: K v (s) — Vj,/[s + 
fj(s)/m]. In other words, the ACF exponentially de- 
cays, but to a non-zero constant, which is the reason 
for nonergodicity and ballistic diffusion. Clearly, this is 
the simplest of possible models, yet physically reason- 
able. The model is non-Markovian, but it allows for a 
three-dimensional Markovian embedding [ItJ : 

x(t) = v(t) 

Q 

mv(t) = -g^V(x) - u(t) - nv(t) + ^2k B Tr 1 ((t) 

u(t) = -vrjv(t) - vu{t) + vy/2k B Trj((t) , (5) 

where £(t) is a zero-centered white Gaussian noise with 
autocorrelation function (C(t)C(t')) = S(t — t'). By inte- 
grating out the auxiliary variable u(t) in the above equa- 
tions (projecting onto the x — v plane), it is not diffi- 
cult to show that this leads to the CLE in Eqs. ((TJ [5]) 
with the memory kernel of Eq. ([3]) , if the initial value of 
u(0) is Gaussian distributed with zero mean and variance 
(w 2 (0)) = k B Trj. This is assumed in the following. The 
initial velocities are thermally distributed. 

We consider the case of a washboard potential with pe- 
riod xo and amplitude strength Vo biased by a constant 
force F, V(x) = ~Vq cos(2ttx / x ) — Fx. It is conve- 
nient to transform Eq. ([5]) into dimensionless quantities 
by scaling time in units of 7 _1 , distance in xq, energy 
in m(a:o7) 2 (which applies for Vq, k B T, and Fxq), u in 
mxo7 2 , and v in 7. We have integrated the system in the 
corresponding nondimensional variables using a standard 
Euler algorithm with time step At = 10~ 4 for v = 0.25, 
Vo = 1, and varying T and F . The behavior of the posi- 
tion variance for 10 particles which started at the origin 
with the velocities thermally distributed is depicted in 
Fig. [T] A striking feature is the regime of intermediate 
hyperdiffusion which ends in the ballistic regime. This 
puzzling behavior, which seems to be rather general, was 
not, however, explained before in physical terms. 

For this we notice that the diffusive behavior can be 
related to the twice-integrated velocity autocorrelation 
function (ACF) (Av(t)Av(t')}. In the stationary limit, 
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FIG. 1: (Color online) Biased diffusion in a washboard po- 
tential with strength Vo = 1 for different values of the bias F 
and the bath temperature T. The mean squared displacement 
exhibits transient hyperdiffusion before it reaches the asymp- 
totically ballistic diffusion regime. The transient hyperdiffu- 
sion is enhanced upon decreasing F and T. The velocities 
are initially thermally distributed, v — 0.25, and n — 10 4 
trajectories are used for the ensemble averaging. 

when the ACF depends on the difference of time argu- 
ments, the normal diffusive behavior emerges when the 
integral of ACF is finite. Its value defines the diffusion co- 
efficient, which is temperature-dependent. For the ther- 
mally distributed velocities and normal diffusion (singu- 
lar limit for v — 0) the ACF decays exponentially to 
zero from (Aw) = ksT/m, with the decay rate 7. This 
yields the Einstein relation, D(T) — k B T/ (rnj). For the 
anomalous GLE, diffusion is described by a spectral bath 
density J{uj) cx j a uj a (0 < a < 2) which corresponds to 
fj(s) = rjas 01 - 1 and yields (Ax 2 (t)) - 2D a (T)t a /T(l + a) 
asymptotically, the integral of the ACF for a ^ 1 ei- 
ther diverges (superdiffusion) , or it tends to zero (ei- 
ther subdiffusion, or bounded motion in trapping po- 
tentials). However, in the absence of forcing, or un- 
der a constant force a generalized Einstein relation al- 
ways holds, D a (T) = k B T / (m^j a ). Furthermore, it is 
easy to show (from the exact expression for K v (s) for 
Vo = 0) that for a > 2, the diffusion is asymptotically al- 
ways ballistic. This is why the occurrence of long-lasting 
hyperdiffusion is rather surprising. For the considered 
model we have (Ax 2 {t)) ~ D^iT^t 2 asymptotically with 
D 2 {T) = k B T/m* , where m* — m(l + 7/1/) is an effec- 
tive mass. This asymptotics holds in the absence of a 
periodic potential, i.e. Vq — 0. 

If one switches on the periodic potential, the ballistic 
diffusion turns over into normal diffusion when F = 
22]. Upon application of a constant force F > 0, the 
particles will start to gradually accelerate when they 
leave the attraction domains of potential wells in the 
phase space due to the random kicks given by ther- 
mal noise. Their initial Maxwellian velocity distribution 
hence will not hold forever. Gradually accelerating, the 
running particles undergo drastic transient heating above 
the bath temperature due to multiple scattering on the 
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FIG. 2: (Color online) The mean square fluctuation of the 
particles velocity for the biased diffusion in Fig. Q] exhibits a 
transient behavior, in which the kinetic temperature Tkm (see 
text) grows from bath temperature to a maximal value T max . 
The transient power law growth, T^ n (t) oc t x , explains the 
hyperdiffusive regime with the corresponding exponent 2 + A 
in Fig. [T] Lower bath temperature and weaker biasing field 
lead surprisingly to stronger heating (see Fig. 



periodic potential. Slower particles, however, are more 
strongly scattered backwards (i.e. are decelerated) by 
the potential wells in comparison with faster ones. This 
in turn yields a growing width (see Fig. [2J of the ve- 
locity distribution, which becomes also skewed towards 
slower particles, see Fig. [31 Indeed, the averaged ki- 
netic energy per particle, K = m(v 2 }/2, can be decom- 
posed as K — K m + Kt- K m = m(v) 2 /2 is the kinetic 
energy associated with mean velocity which asymptoti- 
cally is well described by (v(t)) = Ft/m* when the influ- 
ence of periodic potential becomes negligible. The part 
Kt = m{Av 2 )/2 = /eeTkin/2 is used to define the effec- 
tive kinetic temperature Tid n (t). Importantly, the action 
of F results not only in the growing mean (v(t)), but also 
generates a growing variance of the velocity distribution, 
see in Fig. [2 It is a common practice to characterize 
different sorts of particles with different kinetic temper- 
atures e.g. in plasma physics [28j . In a similar spirit, we 
use a "kinetic temperature" notion which should be used 
with care as it does not correspond to a thermodynamic 
temperature, but rather simply characterizes the width 
of a nonequilibrium velocity distribution, see in Fig. [3j 
Nevertheless, it is a useful concept because it reflects a 
important statistical aspect of the temperature; namely, 
that temperature characterizes the width of the kinetic 
energy distribution. Even if the relative spread of the ve- 
locity distribution around the mean value is rather small 
at the end point of simulations (just a few percents), the 
mean kinetic energy is large and therefore the kinetic 
temperature can also be large, cf. Fig. SJ 

In our numerical experiments, the velocity distribu- 
tion becomes broadened and strongly skewed with the 
"retarded" tail of distribution described by an exponen- 
tial, rather than by a Gauss-Maxwell distribution, cf. 
in Fig. |3J The variance of the distribution grows in 
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FIG. 3: (Color online) Two snapshots of the velocity distribu- 
tion P(v, t), initial and at the end point of simulations in Figs. 
[T][2]for F = 1.5 and T — 0.5. P(v,t) becomes shape invariant 
already for t > 500, after the kinetic temperature reaches its 
maximum, see in Fig. [2] After this happens, only the max- 
imum of the distribution moves accelerating in time. The 
initial distribution is Gaussian with (Aw 2 (0)) = T. The final 
distribution is strongly skewed: its left slope (v < 2999) is 
well fitted by an exponential, while the right slope (v > 2999) 
remains approximately Gaussian. 
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FIG. 4: Dependence of maximal kinetic temperature T max : 
(a) on applied force F for a fixed thermal bath temperature 
T, and (b) on T at a fixed F. 



time in accordance with a power law, T^ n {t) oc t x , cf. 
Fig. [21 This explains the emergence of hyperdiffusion 



(Ax 2 (t)) ~ D 2 (T kin )t 2 = k B T kin t 2 /r 



2+A 



Fig. m 



This regime is, however, only transient. The duration of 
the transient period depends strongly on the potential 
amplitude Vb, the strength of the bias force F, and the 
bath temperature T. 

The hyperdiffusive regime turns over into the ballistic 
diffusion regime when the particles arrive at the maxi- 
mal kinetic temperature T maxi compare Fig. [T]and Fig. 
[2] This nonlinear heating mechanism is quite unusual. 
The heating of plasmas by time- varying stochastic fields 
is well-known. One of the pertinent nonlinear mecha- 
nisms is the so-called Fermi acceleration [2{|. It requires 
but a time-varying driving field (or stochastically oscil- 
lating boundary) . In our case, the "heating" field is how- 
ever constant, and, strikingly enough, the use of weaker 
bias fields heats up the particles ever more strongly, see 
in Fig. [4ja). However, much longer times are required 
then. Moreover, the smaller the bath temperature is, the 
higher is the final kinetic temperature, cf. Fig. 0Jb). 
Both effects are due to the fact that the transient time 
scale becomes longer because the particles take on the 
kinetic energy more slowly in the accelerating field F. 
The corresponding dependencies are stronger than ex- 
ponential. Such a strong sensitivity is surprising. The 
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qualitative physical explanation is as follows: The heat- 
ing is caused by retardation of flying particles when they 
pass over the trapping domains while moving in the bias 
direction. It is appreciably strong as long as the averaged 
energy of the particles does not substantially exceed the 
potential barrier height. For smaller temperatures and 
smaller bias forces it takes an exponentially larger time 
for the particles to escape out of potential wells and to 
arrive on average at a sufficiently high kinetic energy, so 
that their back scattering or deceleration cease to play 
a role. This leads to a larger broadening of the velocity 
distribution [3(j ■ The giant enhancement of ballistic dif- 
fusion reminds one of giant enhancement of normal diffu- 
sion [3ll ] . However, the underlying physical mechanisms 
are quite different. 

The physical systems where the discussed hyperdiffu- 
sive heating effect might be relevant are dusty plasmas 
where heavy tracer particles collide with the gas of light 
particles serving as a thermal bath. Even if our GLE 



description in this case is not directly applicable, there 
exists some partial correspondence between the superdif- 
fusive fractional GLE results and the fractional Kramers 
equation by Barkai and Silbey @. Such a correspon- 
dence is surprisingbecause both descriptions are differ- 
ent, see e.g. in The latter scheme derives from 
linear Boltzmann equations with a fractional scattering 
integral accounting for scattering events which are power- 
law distributed in time. Our physical explanation of the 
transient heating mechanism is more generally applica- 
ble, i.e. it is not restricted by the present GLE model. 
In particular, it is expected to work for the fractional ki- 
netic equations like the fractional Kramers equation, or a 
more general one introduced recently by Friedrich et al. 

We are confident that our work will stimulate fur- 
ther studies and even an experimental validation of this 
intriguing hyperdiffusive behavior. 

This work was supported by the German Excellence 
Initiative via the Nanosystems Initiative Munich (NIM). 
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FIG 1.: Dependence of mean displacement {Ax(t)) on time. 
The parameters are same as in Fig. 1 of main text. 
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FIG 2.: The normalized velocity autocorrelation function as 
a function of time lag r, in the absence and in the presence 
of periodic potential for biased and unbiased diffusion. See 
text for explanation. 

Supplement 

Fig. 1 of this Supplementary Material depicts the de- 
pendence of the mean particle displacement (Ax(t)) on 
time t. First, a regime of transient hyperacceleration 
is clearly seen which accompanies the regime of tran- 
sient hyperdifussion. Second, asymptotically (Ax(t)} is 
well described by (Ax(t)) ~ Ft 2 /(2m*) with the effective 
mass m* — m(l + 7/V) which is m* = 5m for the studied 
parameters. 

Fig. 2 depicts the normalized autocorrelation function 



K v (t,r) = (Av(t+T)Av(t))/(Av 2 (t)) of the velocity fluc- 
tuation, Av(t) = v (t) — (v(t)), as a function of the time 
lag t for t = 500 for three cases: without periodic po- 
tential for a bias F = 1.5 (red line); in the presence of 
periodic potential, Vq = 1, for the same F = 1.5 (blue 
line), and for the same periodic potential in the absence 
of bias, F = 0. In all cases T — 0.5. The analytical 
result in Eq. (4) of the main text (broken line) agrees 
well with numerics (red line). The velocity ACF decays 
exponentially to the constant value v/(y+ r f) = 0.2 which 
indicates nonergodicity, as non-zero correlations hold for- 
ever, for an arbitrary long time lag. Upon switching on 
the periodic potential, the ACF fails to decay (blue line, 
the maximal kinetic temperature T max is already reached 
for studied parameters and t > to = 500). This implies 
that nonergodicity becomes even more pronounced. On 
the contrary, if one switches off the bias F — >• 0, the dif- 
fusion becomes normal and ergodicity is restored in the 
periodic potential. In this case, the ACF decays oscilla- 
tory to zero (green line in Fig. 2). Its integral is finite and 
yields the normal diffusion coefficient (scaled in fcsT/m). 

Next let us discuss shortly the energetics of biased bal- 
listic superdiffusion. After the hyperdiffusion regime is 
over, the particles have arrived at a maximal kinetic tem- 
perature and maximal broadening of the kinetic energy 
distribution. The influence of the periodic potential is 
then negligible. The particles continue, however, to ac- 
celerate in the applied field F, moving as quasi-particles 
with effective mass m* = m(l + 7/^), see in Fig. 1 of 
this Supplementary Material. Asymptotically the mean 
particle position grows like {x(t)) ~ Ft 2 /(2m*) and the 
work spent by F is therefore estimated as W ss F(x) ~ 
F 2 t 2 1 (2m*). Its action is: (i) to accelerate the (real) 
particles to (v) with kinetic energy K m = m(v} 2 /2 — 
(m/m*)W, (ii) to heat them up to K T = m(Av 2 )/2 = 
kBTki n (t)/2. The rest is transferred as an excess heat Q 
to the thermal bath: W = K m + KT + Q. Asymptotically 
Kt -€i W, K m , Q and therefore it can be neglected in 
the energy balance. It thus follows that approximately, 
Q w W - K m w (m* - m)W/m* = jW/{j + v). For 
v 7; most of the acquired work is dissipated as excess 
heat. This is a counter-intuitive conclusion, given that 
the integral friction is zero. For the studied parameters, 
we find that Q « 0.8 W. 



